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^ ■ Abstract 

t— I ■ 

The asymptotic behaviour at late times of inhomogeneous axion-dilaton cosmologies is 
investigated. The space-times considered here admit two abelian space-like Killing vectors. 
These space-times evolve towards an anisotropic universe containing gravitational radiation. 
Furthermore, a peeling-off behaviour of the Weyl tensor and the antisymmetric tensor held 
strength is found. The relation to the pre-big-bang scenario is briefly discussed. 

^ 1 Introduction 

i>: 

The low energy effective action of string theory provides cosmological models that might be 
applicable just below the string scale in the very early universe. In the pre-big-bang scenario 
the universe starts in a low curvature, low coupling regime and then enters a stage of dilaton 

ON ■ driven kinetic inflation [Q] . This picture has been further developed in || . Initially the universe 
consists of a bath of gravitational and dilatonic waves. Some of these collapse leading to the 
^p | birth of a baby inflationary universe. 

The aim here is to investigate the asymptotic future states of a large class of inhomogeneous 
string cosmologies. Late times in the usual general relativistic setting correspond to early times 

. £h ! in the pre-big-bang picture. As will be seen, this provides further evidence that the Milne 
model is not a generic initial state for the pre-big-bang scenario. The issue of a generic initial 
state has also been investigated in || with the emphasis on orthogonal and tilted Bianchi string 
cosmologies. There it was found that a plane wave background is a likely attr actor and not the 
Milne universe as was conjectured in Q. Einstein- Rosen string cosmologies were discussed in 

i §• 

In general relativity the asymptotic evolution of gravitational radiation from a bounded 
source was investigated |7]] and it was found that the Riemann tensor shows a "peeling off" 
behaviour, meaning that different components of the Weyl tensor evolve as different powers of 
the radial coordinate. 

A cosmological peeling off theorem was formulated by Carmeli and Feinstein (8| by using 
earlier results of Stachel || in the case of cylindrical gravitational waves. Cosmological space- 
times of the Einstein- Rosen type can be obtained from these cylindrical gravitational wave space- 
times by interchanging the non-ignorable coordinates. The cosmological peeling off theorem 
decribes the asymptotic future behaviour of inhomogeneous vacuum space-times admitting two 
abelian Killing vectors. The general form of the metric allows space-time homogeneity to be 
broken along one spatial direction, so in general solutions will depend on one spacelike and one 
timelike coordinate. However, this type of metric also encompasses the - by definition - spatially 



homogeneous metrics of Bianchi type I- VII [10] pi. In M it was shown that the Weyl tensor 
shows a peeling off behaviour in terms of the timelike variable. At late times the dominant 
contribution comes from the Weyl scalar ^4 which indicates that the space-time approaches at 
late times a plane wave space-time. 

A similar peeling off behaviour will be found in the case of axion-dilaton string cosmology. 

In the bosonic sector of the low energy effective action of a string theory general relativity 
coupled to two massless scalar fields is obtained in the so-called Einstein frame. 



In four dimensions the action is given in the string frame by [11] 



S = J d^x^-ge-^R + g^d^d^ - -^H a ^H aM ) 



(1) 



where the antisymmetric tensor field strength H a ^ 
Applying the conformal transformation 



dr a Bp^ is introduced. 
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this can be written in the usual Einstein-Hilbert form. 

The physical frame is the string frame since in this frame strings move along geodesies. 
However, since the Weyl tensor is conformally invariant it is justified to use the Einstein frame 
to find solutions and take advantage of the formalism developed for general relativity. 

In the Einstein frame the equations of motion are given by 
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Furthermore in four dimensions the antisymmetric tensor field strength can be written in 
terms of the scalar field, b, as follows 
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2 Equations of motion for axion and dilaton 



Metrics admitting two spacelike abelian Killing vectors can be written in the following general 
form [fill 



2e~ u 
ds 2 = 2e~ M dudv - — j (dx + iZdy) (dx - iZdy) (10) 

M and U are real functions depending on the null variables u and v. Z is a complex function 
depending on u and v. Complex conjugation is denoted by a bar. 

It is convenient to introduce a null tetrad with two real null vectors, 1, n, and two complex 
null vectors which are conjugates of one another, m, m. 



They satisfy the following relations [12] [13] 



l^ = 1 m^Th** = -1 (11) 

and the completeness relation 

9iw = l/j,n u + n^ly - m^fhy - fh^niy. (12) 

Useful formulae for quantities in the Newman-Penrose formalism are given in the appendix. 
Including the axion b and dilaton eft and using Einstein's equations the components of the 
Ricci tensor can be determined. The energy momentum tensor is given by 

T„u = \[M,v + e^b,^ - \g^{^, a ^ + e 2 H a b^)}. (13) 

The following equations of motion are obtained 



2U UU - (U u f + 2M U U U - 4 Z " Z " = <& + e 2 Hl (14) 
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2U VV - (U v ) 2 + 2M V U V - 4 - ^ vZ l„ = (t>l + e 2 H 2 v (15) 
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2Z UV -U U Z V -U V Z U -A-^L = (16) 



2M UV + U U U V - 2 ZuZ " + Z _^v = ^ + e 2<f> buK (17) 
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U uv = U U U V (18) 

2<j) uv -U u (t> v -U v (i) u -2e 24 'b u b v = (19) 

2b uv - U u b v - U v b u + 2(f> u b v + 2(f> v b u = (20) 



Equation (|17| ) is an integrability condition for equations (14) and fllq). The general solution 
of equation (|18|) is given by |1^] 



e- U = f(u)+g(v) (21) 

where / and g are arbitrary functions of their arguments. This can be used to write the 
equations of motion in a more compact way by using / and g as coordinates. 



The equations can be simplified further by introducing a new function S, defined as follows 
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where the prime denotes a derivative w.r.t. to u (v) for / (g). 
So this reduces the system of equations fll4])-(|2"o|) to the following 
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Therefore the components of the Weyl tensor are given by 
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3 Exact solutions 

In this section several exact solutions are discussed. These are of interest since they will already 
show some of the behaviour of the more general space-times which will be discussed in the next 
section. Introducing two new variables, 



f = t- z g = t + z 

and the Ernst potential Z may be conveniently written as 



(31) 
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where £ and oo are real functions. When the imaginary part of the Ernst potential vanishes, 
the two commuting Killing vectors become hypersurface orthogonal. In this case the metric is 
globally diagonalizable. 

In the special case of a diagonal and homogeneous ansatz equation (g|) reads 



C(Ctt + t _1 Ct)-Ct 
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(33) 



which is solved by 



C = af 



(34) 



with a and c constants. 

The equations for the dilaton and axion are given by 



fiu + t-^t-e^bi = 
btt + t~\ + 24> t b t = 



(35) 
(36) 



which is solved by [14| 



b(0 



B 2 
cosh(iVe) + J 1 - -^2 sinh(iV£) 



B cosh(iV£) + Jl-^i sinh(AT£) 



(37) 
(38) 



where dt = td£. This solution can actually be obtained by an SL(2;M) transformation from 
a pure dilaton solution. 

The components of the Weyl tensor are given by 
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Using (|3j) as ansatz for £ and solving (|23j) for cj results in an inhomogeneous solution for u 
]. In this case p^) yields the following two equations 
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For the special value c 



2 this has the solution 
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where O(k) is the step function (0(/c) = for k < 0; 0(«) = 1 for k > 0) and k is an 
arbitrary real parameter. Note that u: can either include F or G, but not both. 

Another solution in this class could be easily generated by applying a mirror transformation. 
Mirror symmetries in the context of metrics with two abelian Killing vectors were discussed 



in 1 15]. The effect of such a transformation in the class of backgrounds considered here is to 
exchange the background metric and the axion-dilaton content. Applying such a transformation 
to the solution above ( pTj ) , ([38]) , (03) , (|4^) results in a solution with a homogeneous dilaton and 
an inhomogeneous, oscillating axion similar to the solution found in ||. Since the solution for 
the axion and dilaton can be obtained by an SL(2;iR) transformation the space-time metric in 
the mirror image is not genuinely diagonal. 

The components of the Weyl tensor are given by 
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For this solution the equations for S, (|2q), (|27|), have the form 
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This system of equations can be solved by introducing a new function £l(f,g) such that 

n f = [u f ] 2 n g = [u g ] 2 (51) 

which implies, using that ojf g = 0, 

n f9 = o. 



This equation has the general solution 

n(f,g) = n 1 (f) + n 2 ( g ) 

and hence using (51) 



Thus S is given by 
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where k is an arbitrary constant. 

With equation (22) the factor e in the expressions for the components of the Weyl tensor 
is given by 



e M = ktf'g'YHl-^e-^. 



(54) 



For 9(«) = 1 and / 
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£ — z fixed, then the components of the Weyl tensor behave as 
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and g = i + z fixed in the limit £ — ► oo. In this case 

the behaviour of ^o an d ^4 is exchanged. In both cases a plane wave space-time is approached. 

However, the plane waves are in opposing directions. 

Furthermore, independent of 0(k), the pre- factor e is given by O {t&~ 3 ) 

Finally, one more class of exact solutions will be discussed. The metric is again assumed to 
-v 
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be diagonal. With C, = e the Ernst equation (|23[) implies 
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This is solved by 
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where Li{} and L 2 {} are arbitrary combinations of the terms in curly brackets. Jo and Yq 
are Bessel functions of the first and second kind. A^, B w , a^ and /3 U are arbitrary constants. 

The first term on its own describes a Kasner solution. The terms containing the Bessel 
functions will be investigated in the next section. The last term introduces the so-called gravi- 
tosolitonic contribution. 

In the simplest case, the gravitosolitonic contribution contains just one term and the constant 
c is set to zero. The constant pre-factor is d. 

Neglecting the contributions from the Bessel functions solutions by Wainwright et al. are 
recovered |16|. These are related to the Kasner metrics by a soliton transformation JlTj |18|. 



Furthermore, for d 2 = a 2 + 3 the metrics are the Ellis-MacCallum solutions of Bianchi type 



III, V and VI [19] [20] [18|. Coordinates i, z adapted to spatial homogeneity are given by |2C][18] 



t = e~ 2m ' z sinh 2mt z = e~ 2mz cosh 2ml 



Depending on the parameter a the models are of Bianchi type III, V or VI. 

As was pointed out by Carmeli et al. [20], in general these space-times are not only singular 



at t = but also at z — ► ±oo. This behaviour depends on the parameters a and d. Only 
the Bianchi models are singularity free at z — > ±oo. If one is interested in a cosmological 
interpretation this seems to confine solutions containing a gravitosolitonic contribution to those 
which are spatially homogeneous. 

Orthogonal Bianchi VI/j models approach a plane wave space-time at future timelike and 



null infinity. The other two approach flat space-time [21] 



4 Asymptotic behaviour of the Weyl tensor and antisymmetric 
tensor field strength 

The Petrov type of the Weyl tensor describes the number of its distinct principal null directions. 



It is determined by the roots of the following equation [13] 

* 4 & 4 + 4f 3 & 3 + 6* 2 & 2 + 4^i6 + ^o = (57) 

In the case of metrics of type (JlQf) the complex Weyl scalars ^1 and ^3 vanish. Therefore 
there are either four distinct or two double roots. So the space-time is of type I or of type D. 



In order to find the Petrov type of a space-time an algorithm developed in [22] can be applied. 



In the case of the metric considered here, the criterion for the space-time being of type D is |2 



@ 



I A = 27 J z (58) 

where I = ^ ^4 + 3^2 and J — ^2(^0^4 - *!)• 

In general the space-time is of type I. However, in the diagonal solution of the last section 
there exist special values for c such that the Weyl tensor is of type D. These are given by c = ±1, 
c = ±Vl-N 2 and c = ±(1 ± ^JA-N 2 ). 

In cosmology one is interested in the evolution of a space-time with time. Therefore future 
timelike and future null infinity will be discussed here. 

Following U Z will be taken of the form 

Z = te- 2 ^{l + i2 X ) (59) 

where ip and x are functions of t and /. In this case tp and x enter more symmetric in the Ernst 
equation and show the same asymptotic behaviour, which in both cases will be 0(t _ a) for tjj 
and x- This asymptotic behaviour show the dilaton and axion as well. This is not surprising 
since the linear terms in the Ernst equation, the equations for the axion and the dilaton have 
the same structure (cf (p3|)-(p5|)). 

At future timelike infinity (t — > 00, z fixed) the Weyl scalars behave as 
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This shows that the space-time does not approach a Kasner model (cf. equations (|39l)-([4l])). 
This can be compared with the future time inifinty of the Doroshkevich, Zeldovich, Novikov 
(DZN) [gj|] |J model for which the components of the Weyl tensor are given by \&o ~ e t , 
^ 2 ~ e M t~ 2 , ^ 4 ~ e M t _1 and the metric is ds 2 = e 2t (dt 2 - dx 2 ) - t d+1 dy 2 - t x - d dz 2 . So 
these models show a similar behaviour of the Weyl scalars at future time infinity. This suggests 
that the nondiagonal Einstein-Rosen space-times approach an anisotropic model similar to the 
DZN model. However, they certainly do not approach the Milne model for which \&q ~ ^2 ~ 
\&4 ~ 0(t~ 2 ). This can be compared with the asymptotic behaviour of diagonal Einstein-Rosen 
axion-dilaton cosmologies. In Q exact solutions were found which approach a DZN universe 
containing axionic, dilatonic and gravitational waves at late times. A similar behaviour is also 
found in electromagnetic Gowdy universes [24]. 

To investigate the behaviour of the Weyl tensor at future null infinity one can use one null 
coordinate, say, / and the timelike coordinate t. 



In these variables the Ernst equation (23) becomes 



z tf + \z tt + Y t ( z t + z } ) = 2(z + zy 



z t z f + -z 2 



(61) 



The linear part of the real and imaginary parts of (|6l|) yield to two cylindrical wave equations 
for ip and x> respectively. These can be solved by Bessel functions (cf (56)). Their asymptotic 
behaviour in the limit t — > oo, / fixed, imply that both ip and \ behave as a(/)t~2, if grav- 
itosolitonic contributions are neglected. These were discussed for diagonal metrics in the last 
section. They seem to lead to a cosmological model only in the spatially homogeneous case. 

In terms of the coordinates t and / the non-vanishing components of the Weyl tensor are 
given by 
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(62) 
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With Z given by ( |59| ) and the asymptotic behaviour of ip and x the components of the Weyl 
tensor behave as follows 



p M 

* ~ —( g 'f[0(t-*)} (65) 

e M 3 

* 2 ~ -i^fg'Pit-*)] (66) 



Ml t l\2 
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* 4 ~ e M (f'y[0(t-2)} (67) 

So as £ — ► oo the component ^4 dominates and the space-time approaches a plane wave 
space-time at future null infinity. 

This behaviour was already found in the exact non-diagonal solution of the last section. 
There the two cases of an ingoing and outgoing wave had been stated explicitly. The same 
applies here as well. Interchanging the null coordinates implies an interchange of the behaviour 
of ^o an( i ^4- Hence outgoing and ingoing waves are exchanged. 

The equations for dilaton and axion are given by 

^ft + ^tt + Y^f + M-^btibf + ^bt] = (68) 

bft + ^b tt + -[b f + b t ]+b t [(l )f + -(f> t } + (f> t [b f + -b t } = 0. (69) 

Hence, in the limit t — > 00 the behaviour of the dilaton and axion is described by, 

0~£>(H) (70) 

6 ~ 0(H). (71) 

This leads to a peeling-off behaviour of the components of the antisymmetric tensor field 
strength H. This is similar to the Einstein-Maxwell case where a peeling-off behaviour for the 
components of the Maxwell tensor is found in addition to that of the Weyl tensor |25| . 

Equation (0) implies that the non-vanishing components of the antisymmetric tensor field 
strength in the null tetrad are given by 

tf(i)(3)(4) = -i\e 2<t> e^f'(b f + l -b t ) (72) 

tf(2)(3)(4) = i^efg'b t . (73) 
This leads to the asymptotic behaviour, 

ff(2)(3)(4) ~ e 2(t> ef0(t-l) (74) 



#(i)(3)(4) ~ e z<f >e-0(t-2). (75) 
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5 Discussion 

The asymptotic behaviour of a large class of inhomogeneous axion-dilaton cosmologies has been 
investigated. At future time infinity the Weyl scalars approach those of an anisotropic space- 
time. At future null infinity the radiative component of the Weyl scalars dominates. This might 
be interpreted as that this class of inhomogeneous string cosmologies approaches an expanding 
anisotropic background on which gravitational waves propagate. Axion and dilaton become 
negligible at late times. 

A peeling-off behaviour of the components of the Weyl tensor has been found similar to 
that in the case of vacuum general relativity. The components of the antisymmetric tensor field 
strength also show a peeling-off behaviour. This is analog to the asymptotic behaviour of the 
Maxwell tensor in the Einstein-Maxwell case. 

The future asymptotic state in general relativity can be interpreted as the past asymptotic 
state for the pre-big-bang scenario. In the class of inhomgeneous metrics considered here, this 
means that an anisotropic model and not the Milne universe is a likely past attractor. 

Furthermore, the type of metric which was studied here might also be relevant to the recently 
proposed pre-big-bang bubble picture. Together with appropriate boundary conditions metrics 
of the general form (|l(]) also describe colliding plane wave space-times. This might be used to 
find a "dynamical" description of the formation of pre-big-bang bubbles p| . 
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Appendix: Components of the Weyl and Ricci tensors 

Here all the non-vanishing components of the Weyl tensor and the Ricci tensor in the Newman- 
Penrose formalism are given. The signature of the metric is (H — -) and Einstein's equations 
are given by 

Rfj,u — Tid^R = —Tfja,. (76) 

The tetrad metric to raise and lower terad indices is given by 
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Tetrad indices are enclosed in brackets, which will be omitted if it is unambiguous. Tetrad 
indices run from 1 to 4. 

For a space-time with metric ( |To| ) a null tetrad basis is provided by (I! 
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The Newman-Penrose spin coefficients have been calculated using Maple [26], 
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F u denotes the partial derivative *L-. 

In the Newman-Penrose formalism the components of the Weyl tensor are given by five 



complex scalars [13| 
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and the components of the Ricci tensor are denoted by the following four real and three 



complex scalars [13| 
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Using the relations given in [ 13 1 the components of the Weyl and the Ricci tensor are derived. 
The components of the Weyl tensor are given by 



M 



$0 = 


e 


(Z + Z) 2 


^2 = 


- -r M ( 2C 


^4 = 


e M 


(Z + Z) 2 


$1 = 


= ^3 = 



(Z + Z)(Z VV + Af w Z„ - U V Z V ) - 2{Z V 



Z U Z V 



2A 



(Z + Z)(Z UU + M U Z U - £/ u Z u ) - 2(Z U ) 2 



and the non- vanishing components of the Ricci tensor are 



(82) 
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$00 
$02 
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2U VV - (U v ) 2 + 2M V U V - 4 
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(Z + Zf 



1 e 



7V7 



2(Z + Z) 2 



[(Z + Z)(2Z TO - c^z, - tr„z u ) - 4Z U Z 1) ] 



2M utI + £7^ - 2 



Z U Z V + Z U Z V 

~JzVzy~ 



e M 

$22 = -r 



2C/„„ - (C4) 2 + 2[/ u M„ - 4 Z "''" 



A 



" 24 



2M UV + 4E/ U „ - 3U U U V - 2 



(z + zy 

Z U Z V + Z U Z V 



(83) 



(Z + Z) 2 
The non- vanishing components and the trace in a space-time with metric ( |l0|) are given by 

Too 



M-U 



T- 



Z-Z, 



22 



Y^zw»b> + ^^ T2 3 = 2 eM_t/ zTz 



[33 = e 



M-U 



zz 

z + z^ 



u (f) v + e ^6,;] 
^ + e 2 *6 u 6„] T = -e M [<t> u <f> v + e^b u b v } (84) 
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and using Einstein's equations the only non-vanishing components of the Ricci tensor are as 
follows 

Roo = -\\<& + ^ h l\ R oi = -5&A + e 2 H u b v ] R u = -\[4>l + e 2 Hl\. (85) 

So the non- vanishing Newman- Penrose scalars <& a b and A are given by 

M 

$00 = -j- [<f% + e 2 %] $11 = V [^ + e 2* bu6t; ] 

$22 = — [<t>l + e 2 H 2 u ] A = -%z[K4>v + e^bub,,]. (86) 
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